Blocking the gate to ligand entry in human hemoglobin

His(E7) to Trp replacements in HbA lead to markedly biphasic bimolecular CO rebinding after laser photolysis. For isolated mutant subunits, the fraction of fast phase increases with increasing [CO], suggesting a competition between binding to an open conformation with an empty E7 channel and relaxation to blocked or closed, slowly reacting states. The rate of conformational relaxation of the open state is ～18,000 s(-1) in α subunits and ～10-fold faster in β subunits, ～175,000 s(-1). Crystal structures were determined for tetrameric α(WT)β(Trp-63) HbCO, α(Trp-58)β(WT) deoxyHb, and Trp-64 deoxy- and CO-Mb as controls. In Trp-63(E7) βCO, the indole side chain is located in the solvent interface, blocking entry into the E7 channel. Similar blocked Trp-64(E7) conformations are observed in the mutant Mb crystal structures. In Trp-58(E7) deoxy-α subunits, the indole side chain fills both the channel and the distal pocket, forming a completely closed state. The bimolecular rate constant for CO binding, k'(CO), to the open conformations of both mutant Hb subunits is ～80-90 μm(-1) s(-1), whereas k'(CO) for the completely closed states is 1000-fold slower, ～0.08 μm(-1) s(-1). A transient intermediate with k'(CO) ≈ 0.7 μm(-1) s(-1) is observed after photolysis of Trp-63(E7) βCO subunits and indicates that the indole ring blocks the entrance to the E7 channel, as observed in the crystal structures of Trp(E7) deoxyMb and βCO subunits. Thus, either blocking or completely filling the E7 channel dramatically slows bimolecular binding, providing strong evidence that the E7 channel is the major pathway (≥90%) for ligand entry in human hemoglobin.     

Local selection across a latitudinal gradient shapes nucleotide diversity in balsam poplar, Populus balsamifera L

Molecular studies of adaptive evolution often focus on detecting selective sweeps driven by positive selection on a species-wide scale; however, much adaptation is local, particularly of ecologically important traits. Here, we look for evidence of range-wide and local adaptation at candidate genes for adaptive phenology in balsam poplar, Populus balsamifera, a widespread forest tree whose range extends across environmental gradients of photoperiod and growing season length. We examined nucleotide diversity of 27 poplar homologs of the flowering-time network-a group of genes that control plant developmental phenology through interactions with environmental cues such as photoperiod and temperature. Only one gene, ZTL2, showed evidence of reduced diversity and an excess of fixed replacement sites, consistent with a species-wide selective sweep. Two other genes, LFY and FRI, harbored high levels of nucleotide diversity and exhibited elevated differentiation between northern and southern accessions, suggesting local adaptation along a latitudinal gradient. Interestingly, FRI has also been identified as a target of local selection between northern and southern accessions of Arabidopsis thaliana, indicating that this gene may be commonly involved in ecological adaptation in distantly related species. Our findings suggest an important role for local selection shaping molecular diversity and reveal limitations of inferring molecular adaptation from analyses designed only to detect species-wide selective sweeps.     

Complete genome sequence of Methanosaeta concilii, a specialist in aceticlastic methanogenesis

The genome sequence of the aceticlastic methanoarchaeon Methanosaeta concilii GP6, comprised of a 3,008,626-bp chromosome and an 18,019-bp episome, has been determined and exhibits considerable differences in gene content from that of Methanosaeta thermophila.     

The Neovolcanic Axis Is a Barrier to Gene Flow among Aedes aegypti Populations in Mexico That Differ in Vector Competence for Dengue 2 Virus

Background

Aedes aegypti is the main mosquito vector of the four serotypes of dengue virus (DENV). Previous population genetic and vector competence studies have demonstrated substantial genetic structure and major differences in the ability to transmit dengue viruses in Ae. aegypti populations in Mexico.

Methodology/Principal Findings

Population genetic studies revealed that the intersection of the Neovolcanic axis (NVA) with the Gulf of Mexico coast in the state of Veracruz acts as a discrete barrier to gene flow among Ae. aegypti populations north and south of the NVA. The mosquito populations north and south of the NVA also differed in their vector competence (VC) for dengue serotype 2 virus (DENV2). The average VC rate for Ae. aegypti mosquitoes from populations from north of the NVA was 0.55; in contrast the average VC rate for mosquitoes from populations from south of the NVA was 0.20. Most of this variation was attributable to a midgut infection and escape barriers. In Ae. aegypti north of the NVA 21.5% failed to develop midgut infections and 30.3% of those with an infected midgut failed to develop a disseminated infection. In contrast, south of the NVA 45.2% failed to develop midgut infections and 62.8% of those with an infected midgut failed to develop a disseminated infection.

Conclusions

Barriers to gene flow in vector populations may also impact the frequency of genes that condition continuous and epidemiologically relevant traits such as vector competence. Further studies are warranted to determine why the NVA is a barrier to gene flow and to determine whether the differences in vector competence seen north and south of the NVA are stable and epidemiologically significant.     

Shape complementarity of protein-protein complexes at multiple resolutions

Biological complexes typically exhibit intermolecular interfaces of high shape complementarity. Many computational docking approaches use this surface complementarity as a guide in the search for predicting the structures of protein-protein complexes. Proteins often undergo conformational changes to create a highly complementary interface when associating. These conformational changes are a major cause of failure for automated docking procedures when predicting binding modes between proteins using their unbound conformations. Low resolution surfaces in which high frequency geometric details are omitted have been used to address this problem. These smoothed, or blurred, surfaces are expected to minimize the differences between free and bound structures, especially those that are due to side chain conformations or small backbone deviations. Despite the fact that this approach has been used in many docking protocols, there has yet to be a systematic study of the effects of such surface smoothing on the shape complementarity of the resulting interfaces. Here we investigate this question by computing shape complementarity of a set of 66 protein-protein complexes represented by multiresolution blurred surfaces. Complexed and unbound structures are available for these protein-protein complexes. They are a subset of complexes from a nonredundant docking benchmark selected for rigidity (i.e. the proteins undergo limited conformational changes between their bound and unbound states). In this work, we construct the surfaces by isocontouring a density map obtained by accumulating the densities of Gaussian functions placed at all atom centers of the molecule. The smoothness or resolution is specified by a Gaussian fall-off coefficient, termed "blobbyness." Shape complementarity is quantified using a histogram of the shortest distances between two proteins' surface mesh vertices for both the crystallographic complexes and the complexes built using the protein structures in their unbound conformation. The histograms calculated for the bound complex structures demonstrate that medium resolution smoothing (blobbyness = -0.9) can reproduce about 88% of the shape complementarity of atomic resolution surfaces. Complexes formed from the free component structures show a partial loss of shape complementarity (more overlaps and gaps) with the atomic resolution surfaces. For surfaces smoothed to low resolution (blobbyness = -0.3), we find more consistency of shape complementarity between the complexed and free cases. To further reduce bad contacts without significantly impacting the good contacts we introduce another blurred surface, in which the Gaussian densities of flexible atoms are reduced. From these results we discuss the use of shape complementarity in protein-protein docking.     

Purification and multimer formation of allurin, a sperm chemoattractant from Xenopus laevis egg jelly

Allurin, a sperm chemoattractant isolated from Xenopus laevis egg jelly, can be purified in one step from an extract of diffusible jelly proteins ("egg water") using a FPLC or HPLC anion exchange column and a multi-step NaCl gradient. Allurin homomultimers were detected by Western blotting with antibodies prepared against the purified protein or peptides within the protein. Allurin multimers were stable and resisted dissociation by SDS and beta-mercaptoethanol. Alkylation of allurin provided evidence for two free sulfhydryl groups but did not eliminate multimer formation, suggesting that intermolecular disulfide bond formation is not required for allurin aggregation. Concentration of egg water was accompanied by a reduction of chemoattractant activity that could not be fully accounted for by homomultimer formation. Rather, the presence of a multiphasic dose-activity curve upon partial purification and formation of hetero-allurin complexes during concentration suggested that egg water may contain allurin-binding proteins that reduce multimer formation and activity.     

The tick saliva immunosuppressor, Salp15, contributes to Th17-induced pathology during Experimental Autoimmune Encephalomyelitis

Salp15 is a tick saliva protein that inhibits CD4⁺ T cell differentiation through its interaction with CD4. The protein inhibits early signaling events during T cell activation and IL-2 production. Because murine Experimental Autoimmune Encephalomyelitis development is mediated by central nervous system-infiltrating CD4⁺ T cells that are specific for myelin-associated proteins, we sought to determine whether the treatment of mice with Salp15 during EAE induction would prevent the generation of proinflammatory T cell responses and the development of the disease. Surprisingly, Salp15-treated mice developed more severe EAE than control animals. The treatment of EAE-induced mice with the tick saliva protein did not result in increased infiltration of T cells to the central nervous system, indicating that Salp15 had not affected the permeability of the blood-brain barrier. Salp15 treatment did not affect the development of antibody responses against the eliciting peptide or the presence of IFNγ in the sera. The treatment with Salp15 resulted, however, in the increased differentiation of Th17 cells in vivo, as evidenced by higher IL-17 production from PLP_139-151-specific CD4⁺ T cells isolated from the central nervous system and the periphery. In vitro, Salp15 was able to induce the differentiation of Th17 cells in the presence of IL-6 and the absence of TGFβ These results suggest that a conductive milieu for the differentiation of Th17 cells can be achieved by restriction of the production of IL-2 during T cell differentiation, a role that may be performed by TGFβ and other immunosuppressive agents.     

Competing RNA secondary structures are required for mutually exclusive splicing of the Dscam exon 6 cluster

Alternative splicing of eukaryotic pre-mRNAs is an important mechanism for generating proteome diversity and regulating gene expression. The Drosophila melanogaster Down Syndrome Cell Adhesion Molecule (Dscam) gene is an extreme example of mutually exclusive splicing. Dscam contains 95 alternatively spliced exons that potentially encode 38,016 distinct mRNA and protein isoforms. We previously identified two sets of conserved sequence elements, the docking site and selector sequences in the Dscam exon 6 cluster, which contains 48 mutually exclusive exons. These elements were proposed to engage in competing RNA secondary structures required for mutually exclusive splicing, though this model has not yet been experimentally tested. Here we describe a new system that allowed us to demonstrate that the docking site and selector sequences are indeed required for exon 6 mutually exclusive splicing and that the strength of these RNA structures determines the frequency of exon 6 inclusion. We also show that the function of the docking site has been conserved for ~500 million years of evolution. This work demonstrates that conserved intronic sequences play a functional role in mutually exclusive splicing of the Dscam exon 6 cluster.     

Introduction of a disulfide bond leads to stabilization and crystallization of a ricin immunogen

RTA1-33/44-198 is a catalytically inactive, single-domain derivative of the ricin toxin A-chain (RTA) engineered to serve as a stable protein scaffold for presentation of native immunogenic epitopes (Olson et al., Protein Eng Des Sel 2004;17:391-397). To improve the stability and solubility of RTA1-33/44-198 further, we have undertaken the design challenge of introducing a disulfide (SS) bond. Nine pairs of residues were selected for placement of the SS-bond based on molecular dynamics simulation studies of the modeled single-domain chain. Disulfide formation at either of two positions (R48C/T77C or V49C/E99C) involving a specific surface loop (44-55) increased the protein melting temperature by ~5°C compared with RTA1-33/44-198 and by ~13°C compared with RTA. Prolonged stability studies of the R48C/T77C variant (> 60 days at 37°C, pH 7.4) confirmed a > 40% reduction in self-aggregation compared with RTA1-33/44-198 lacking the SS-bond. The R48C/T77C variant retained affinity for anti-RTA antibodies capable of neutralizing ricin toxin, including a monoclonal that recognizes a human B-cell epitope. Introduction of either R48C/T77C or V49C/E99C promoted crystallization of RTA1-33/44-198, and the X-ray structures of the variants were solved to 2.3 A or 2.1 A resolution, respectively. The structures confirm formation of an intramolecular SS-bond, and reveal a single-domain fold that is significantly reduced in volume compared with RTA. Loop 44 to 55 is partly disordered as predicted by simulations, and is positioned to form self-self interactions between symmetry-related molecules. We discuss the importance of RTA loop 34 to 55 as a nucleus for unfolding and aggregation, and draw conclusions for ongoing structure-based minimalist design of RTA-based immunogens.